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We examine quantum normal typicality and ergodicity properties for quantum systems whose dynamics are 
generated by Hamiltonians which have residual degeneracy in their spectrum and resonance in their energy gaps. 
Such systems can be considered atypical in the sense that degeneracy, which is usually a sign of symmetry, 
is naturally broken in typical systems due to stochastic perturbations. In particular, we prove a version of 
von Neumann's quantum ergodic theorem, where a modified condition needs to hold in order to have normal 
typicality and ergodicity. As a result, we show that degeneracy of spectrum does not considerably modify the 
condition of the theorem, whereas the existence of resonance is more dominant for obstructing ergodicity. 
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I. INTRODUCTION 


Statistical mechanics has proved to be a successful theory for macroscopic systems. One of the cornerstones of statistical 
mechanics is the ergodic hypothesis, which colloquially states that the fraction of time state of a systems spends in a given 
subspace of its state space is proportional to the fraction of the surface occupied by this subspace [1-4]. The situation in 
quantum statistical mechanics becomes specially interesting in light of basic differences with classical statistical mechanics. 
Such differences are caused by the very mathematical structure of quantum mechanics, which allow “typicality” behaviors to 
emerge [5-10], 

The first, seminal attempt to put the ergodic hypothesis in the form of a rigorous theorem was made by von Neumann through 
proving a “quantum ergodic theorem” [11, 12], The validity of this theorem, however, was heavily debated later in the literature 
[13-15], It was just recently that this theorem was revisited carefully in Ref. [16] where the earlier criticisms were refuted, 
and it was rigorously proved that the theorem (with the set of sufficient conditions assumed by von Neumann) is indeed valid. 
Specifically, it has been shown that von Neumann’s statement is in fact a more general property than the ergodicity, called normal 
typicality. This reaffirmation of the quantum ergodic theorem resolves a long-standing issue with quantum statistical mechanics, 
and is of fundamental importance. Interest in studying underlying laws of thermodynamics and statistical mechanics has been 
recently reinvigorated (mostly by advances in quantum information science), which also includes revisiting (emergent) properties 
such as ergodicity, thermal equilibration, or out-of-equilibrium fluctuations [17-34]. Thus a clear form of the quantum ergodic 
theorem can enrich such attempts to better understand the principles of statistical mechanics in light of quantum mechanics. 

Here we further the proof of Ref. [16], and lift some necessary conditions of the quantum ergodic theorem, and investigate 
how this may affect the theorem. Non-degeneracy of energy spectrum and energy gaps both have been assumed in the proof 
of the quantum ergodic theorem. Non-degeneracy and non-resonance are mainly due to small and uncontrollable interactions 
of quantum systems with their environment, and thus exist in typical systems. Here we, however, consider atypical quantum 
systems whose dynamics are generated by Hamiltonians which have residual degeneracy in their spectrum (not completely lifted 
by external perturbations) and resonance in their energy gaps. We find conditions under which normal typicality and ergodicity 
may still hold to some extent. 

The structure of this paper is as follows. In Sec. II, we set notations and recall necessary definitions and results from literature. 
In Sec. Ill, we state our main result in Theorem 1, which shows how non-degeneracy and non-resonance conditions modify 
the condition for normal typicality and ergodicity. We next provide the proof of Theorem 1 in Sec. IV. The paper ends by a 
discussion and summary of our result, followed by two appendices where some necessary details are proved. 


II. PRELIMINARIES 


In this section, we set up some preliminaries and remind pertinent definitions and results (see Ref. [16] for details). 
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FIG. 1. Schematic of the energy, gap, and energy sum structures. Each energy eigenvalue E a has degeneracy e Q (a £ {1,2,. .. , De}); 
each gap value G k = E a — Ej (for some k(a,7)), where a 7 ^ 7, has degeneracy and each energy sum value F m = E a + E 7 
(for some m(a,7)) has degeneracy / m . Here Q 0 denotes the vanishing-gap set, which has g 0 = He elements. In addition, we have 
9 k = De{De — 1). We also define Dg = maxk^o 3 k and Dp = max m f m . The case of the non-resonant Hamiltonian is given by 
Dg = 1, or equivalently, Df = 2 (see Appendix A). 


A. Hamiltonian and dynamics 

Let us assume that we have a quantum system whose associated Hilbert space is 9~( (where D = dim(.?/")). The dynamics of 
this system is generated by the Hamiltonian 


d b ^d 

H= J 2 d) 

c*=1 

where E a is an e Q -fold degenerate energy eigenvalue (E a 7 ^ Ep iff a / ji) and n a = \a, a)(ct,a| represents the 

corresponding eigenprojection. Here e « = E>, in which Dp is the number of distinct eigenvalues. Any initial state p(0)} 

evolves into 


^(t)) = e ITH \ip( 0)), (2) 

after time t, where we have assumed h=l (here and hereafter). 

We also need to define (unnormalized) vectors | <p a ) and (normalized) vectors \ip a ) belonging to the nth energy shell as 

e a 

n a |v>(o)) = y^5 aa |Q, a) 

a—1 

= \v a ) = C a \ip a ), ( 3 ) 

where s aa = (a, a\i/>(0)} and c a is the normalization factor (|c a | 2 = (£> a | 'fa})- 


B. Gap and energy-sum structures 

We set the shorthand Eg — E a =: Gk to denote the energy gap, for some k. For a given value of Gk, there may exist several 
ordered pairs (a,/3) for which Ep — E a equals this given value of G. Thus we define := {(a,/3) | Ep — E a = Gk}- It 
is evident that \Qk\ = ,9k- where is the degeneracy of the energy gap Gk, and A denotes cardinality of set A. Similarly, 
for sum F m of distinct energies, we define := {( 01 , 7 ) | E a + E^ = H m } and / m = |{F m |. For later use, we also define 
Dq = maxk^o 9k and Dp = max m / m . Figure 1 illustrates a schematic of the gap and energy-sum structures. 

A Hamiltonian is called non-degenerate if none of its energy eigenvalues E a is degenerate, i.e., e a = 1 Va, or equivalently 
Dp = D. Similarly, a Hamiltonian is called non-resonant if none of its energy gaps is degenerate, i.e., 9 k = 1 Vk 7 ^ 0, or 
equivalently Dq = 1- Later in Appendix A we show that the non-resonance condition implies f m = 2 Vm, or equivalently Dp = 
2. Thus, in general we have Dp A 2. We were, however, not able to find a general relation between the gap-related quantity 
Dq and the corresponding energy-sum-related quantity Dp. But intuitively it seems that a relatively small Dq corresponds to a 
relatively small Dp. 
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FIG. 2. Schematic structure of the Hilbert space 9~l = decomposed to M “macrospaces” or “Gibbs cells” (green cells) according 

to a complete set of macro-observables Each (experimentally accessible) Gibbs cell includes many “microstates” (gray cells), as 

implied by the assumption dim(d~C) = <C dim (9i) = D, Mv. 


C. Measurements, macrostates, and microstates 

We partition 9~( according to a complete set of given (but otherwise arbitrary) orthogonal projections { P v }^'Li , in which P v 
has rank d v and -P^ = 1 [ 11], as 

rt = (4) 

with = {P v \v), |tt) £ 9{}. It is evident that 

dim(yC) = d v , (5) 

M 

J2d, = D. ( 6 ) 

l/=l 

The set {P V }^L 1 can be associated with a complete measurement on the system [35]; or equivalently, each can be 
interpreted to represent a “macrospace/macrostate” or “Gibbs cell” of the system in which all “microstates” yield the same result 
for the measurement—see Fig. 2. 


D. Averaging and von Neumann’s lemmas 

Let E\-[/(x)] = J 2) /( x) px(x)dx denote the average or expected value of f{x) when A is a random variable defined on a 
set T> C ]R with a given probability distribution p x {x). In addition, the variance of A' is given by Vx [x] = Ex [(a-’ ~ Ex [x]) 2 ] 
[3]. From positivity of the variance, Vx[i] 0, we immediately have 

ExN ^ Ex[x 2 ]. (7) 

Similarly, from Ex[(a?i — a^) 2 ] ^ 0 for two random variables X\ and x-i, we have 

2¥, x [x\X2] < E.y[ii] + ¥. x [x\). (8) 

Hence, when X\ and x 2 have identical expected values and variances, we have 

Ea [x\x 2 ] < E^[ati] + Y x [x{\ . (9) 

One can also see from the definition of the variance that it satisfies the following property: 

K K K 

Vx[5>i 

i— 1 * = 1 

where C x [x{, Xj] = Ex[ij] Ex[ij] — E\'[xjXj] is the covariance of the random variables x t and x :) [2, 3], 

In addition to the basic definitions above, we shall also need a result known as Markov’s inequality [3], 

Prob[X ^ B] < E x [x\/B, B > 0. (11) 

From this we say X ^ B holds for (1 — <5) most A's if ProbfAT ^ B] ^ 1 — <5 [16]. 

The following lemmas—proven by von Neumann (adopted from Ref. [16])— will also be essential later: 


= y^Vx[:ri]+ ^ Cx[a 


2 5 5 


( 10 ) 
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Lemma 1 There exists a number C > 1 such that when two natural numbers d and D satisfy the condition 

\ogD_ d 1 

D D C ’ 


/or uniformly-randomly distributed unitaries U £ SU(D) (according to the Haar measure [2, 36]) we have 


E r 


max 
U 7 ye{i,...,D} 


E U ^ U > 


(j 


t=i 


2 t ^ log-D ^ 1 
< D < ^2’ 


E r . 


max 




d \ 2 1 9d log D 


l=\ 




D 2 


< 


C 3 ' 


( 12 ) 


(13) 

(14) 


Lemma 2 Let |£) be a uniformly-randomly distributed state from a D-dimensional Hilbert space and P be a rank-d projection 
thereon. Then 


E?[i™il 2 ] 

v £ [||P|OII 2 ] 


d_ 

D 1 

1 f d\ 2 (D-d) 
d\Dj (£> + !)’ 


(15) 

(16) 


For the sketch of proof see Appendix B [2], 


E. Long-time average and long run 

For a dynamical function Y (r), we define the long-time average Y as follows [19]: 

Y = lim * [ Y (r)dr. (17) 

T —>-oo T J o 

We also say that a statement S(t) holds for (1 — S') fraction of the time in the long run when 

lim inf ^|{r G [0,Tl; S(t) holds}! ^ (1 — S'), (18) 

where | A\ indicates the Lebesgue measure or size of set A C R [16]. 

Remark. —In the standard context of the quantum ergodic theorem we focus on the T —> oo time averages. It is, however, 
interesting to see how the setting of the theorem and associated results may change for finite-time averages. This subject is 
beyond the scope of the current paper. See, for example. Refs. [23, 24, 30, 33, 34] for discussions on finite time scales for 
equilibration. 


F. Normal typicality 

The normal typicality property in the sense of von Neumann implies that for almost all partitionings of a Hilbert space 
Of according to a complete set of M rank-<i„ orthogonal projections (say { P v }^Lf) as in Eq. (4), and for all initial states 
|t/>(0)) £ Of, a system is e-6' normal if for (1 — S') fraction of the time in the long run we have [16] 

\\\P^( T ))\f-d v /D\^-^^, Mv £ (19) 

where ||iz|| = sj (w|tt)' (for |u) £ Of) indicates the Euclidean vector norm. 

A sufficient condition for this result to hold is that 

L = |||P,|/-(r))|| 2 - d v /D\ 2 < 6'(^) 2d £, Mv £ (20) 

in the long run [16]. The proof is immediate by contradiction. To see this, assume Y (r), p, and q to be the left-hand side (LHS) 
of Eq. (19), EyJa v /(MD) , and S'/M, respectively. Let us assume that Y(t) ^ p (^ 0) for at least q fraction of the total time 
T; this assumption is clearly a violation of Eq. (19). Hence the time average of Y 2 (t) on the interval [0, T] is / qp 2 ; which in 
turn violates Eq. (20). Thus Eq. (20) is a sufficient condition for the property (19). 
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G. (Relative) Ergodicity 


We call a closed quantum system evolving via a time-independent Hamiltonian e-ergodic relative to a given set of complete 
observables (measurements) {P V }^L 1 if for every initial state |-0(O)) we have 


IIW ( t )>|| 2 



W G {1,..., M}. 


( 21 ) 


In other words, when the RHS of Eq. (21) becomes negligible, we conclude that the average time the system spends in a 
certain macrostate (induced by the P v measurements) is proportional to the relative size of that Gibbs cell (here denoted by 
d v ) compared to the size of the entire available Hilbert space (here denoted by D). If the above property holds for (almost) all 
possible measurements, we call the system ergodic. This statement somehow resembles the traditional statement of the ergodic 
hypothesis in classical statistical mechanics [3], 

There is a clear connection between normal typicality and ergodicity. We note that 


\P,mr)W - d v /D\ = ||P,|V'(r))|| 2 - d v /D\ 


= lim —- 

T—foo T 2 


\ P v\'f( T ))\\ 2 - d„/D) dr 


Eq. (7); -Y ^.Y 2 l 

^ lim — 

T—foo T 


= \\\PM{TW~d„/DY 


\ p v\4 , (. T )}\\ 2 - d v /D\ dr 


or equivalently. 


\P^{.rW-d v /D \ 2 ^L. 


( 22 ) 

(23) 


Hence normal typicality (i.e., a relatively small upper bound on L) implies ergodicity [16]. We shall use this fact later. 


III. MAIN RESULT 


Our aim is to see whether and how the results of Refs. [12, 16] are modified in the case of degenerate and resonant Hamil¬ 
tonians. The following theorem, a modified version of von Neumann’s quantum ergodic theorem, encapsulates our main result. 
Note that a stronger (or perhaps tighter) version of the theorem may still be possible, but we shall not discuss it here. 

Theorem 1 Let 0 < 6. S' f 1 and e > 0, and consider a D-dimensional Hilbert space Of associated with a quantum system. 
This system is £-8' normal for (1-8) fraction of all Hilbert-space decompositions Of = ®^f =1 0f u , with given (fixed) dimensions 
{d u } induced by measurements, if 


max 


r iom 2 

\ Cl 88'e 2 


( (D F -2)dl \ 

v 10.D log D ) 


log D ^ d u ^ 1 
D < If < C’ 




(24) 


where C is a number satisfying Lemma 1, and Dp is the maximum degeneracy of the energy sum structure of the system (defined 
in Subsec. IIB). 


In the absence of resonance (the case where Dp = 2), Eq. (24) implies that as long as the dimension of each Gibbs cell 
is sufficiently large (d v 1), yet sufficiently small compared to the dimension of the entire Hilbert space (d v <C D), the 
system will feature normal typicality for a certain fraction of the measurements. This is the standard quantum ergodic theorem 
[16]. Note, however, that in Eq. (24) degeneracy does not show up explicitly, whereas resonance appears to be important for 
obstructing or allowing normal typicality (and ergodicity)—because of the existence of the factor (Dp — 2)d 2 /(10Z?logD). 
From this one can conclude that the impact of resonance would be relatively small if 

Dp < (101ogD/D)M 2 , (25) 

where we have used the crude estimate d v « 0(D/M ) [see Eq. (6)]. 

We also remark that for a given pair (d„, D) there may exist several G s satisfying Lemma 1. In this case we choose the largest 
one. As an example, for a system of 10 2 spin-1/2 particle with D = 2 10 « 10 30 , when d u « 10 8 , and M « 10 22 , we have 
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log D/D ft 10 29 , dv/D ft 10 22 . Thus, we can choose C ft 10'. Additionally, one may be able to choose e and such that 
they satisfy the condition (24), for sufficiently negligible values of 8, and hence observe ergodicity according to Eqs. (20) and 
( 21 ). 

It would be instructive to consider specific physical examples (e.g., spin chains with strong symmetries and too large degenera¬ 
cies and resonance) in which one or parts of the conditions of Theorem 1 do not hold, whereas one could explicitly demonstrate 
that the system is not ergodic. We, however, leave this (important) task for future investigations, and here focus mainly on the 
formalism of the quantum ergodic theorem. 


IV. PROOF OF THEOREM 1 


The logic of the proof goes as follows. We find an upper bound X on L. In order for normal typicality (and whence— 
according to Eq. (23)— ergodicity) to hold, we shall require that this upperbound A' itself be not greater than the value given by 
the sufficient condition for e-8' normal typicality (20), specifically, the RHS of this equation: B = 5 '(e/ M) 2 (d v / D). That is, 

Eq. (23) 

ergodicity [LHS of Eq. (21)] < L^X^B [RHS of Eq. (20)]. (26) 

Next we shall employ Markov’s inequality and require that its RHS be in turn upperbounded by 8, 

Prob[X > B\ < < 8 , (27) 

jD 

which is tantamount to 

Prob[A ^ B] Eq ' = i 26 ^ p r ob[ e-8' normal typicality ] = Prob[ ergodicity ] > 1 — 8. (28) 

This implies that e-8' normal typicality and ergodicity hold for (1 — 8) most measurements { P v } with fixed ranks { d v }. We note 

that in order for this construction to work, we need to consider how our parameter d v should be appropriately chosen. There 
are two relevant conditions; the first one comes from the condition on d„ to apply Lemma 1, and the second condition comes 
from the requirement in the last part of Eq. (27). We combine these conditions to find how/when an e-8' (with given values) 
normality/ergodicity is implied. In the sequel, we follow the steps elaborated here. 

Expanding L [Eq. (20)] in terms of {] 8p a )} yields 

L ={d v /Df - 2 (djD) (8p a \P u \v p ) + ]T e~ iT ( E P- E °‘ +E 'i- E ^{(p a \P u \tpp). (29) 

a/3 a/37 <t 

In the following, we consider this expression term by term as 

L = {d v /D) 2 + L 2 +L 3 . (30) 


The second term (L 2 ) concerns degeneracy. If the Hamiltonian is non-degenerate, L 2 can be simplified by using e lT{ E “ E,y > 

8 a(H tO 


D 

X ( 2 nd) = -2(d„/Z>) p u\!p a ). (31) 

a—1 

In the presence of degeneracy, we have 

£ (d) = —2(d v /D) e.~ iT ^ E f>~ Eo ‘\8p a \P v \ipp) 
a/3 

Eq = (3) -2(d„/Z>) ^ 8 EaEf ,s* aa sp b (a,a\P„\p,b) 

a(3,ab 

= - 2 {d v /D) ^ s* aa s ab (a,a\P„\a,b) 

cx,ab 
D e 

= - 2 (d v /D) ^2 \c a \ 2 (<p a \P„\<p a ) 

a—1 

D e 

= - 2 (d v /D) ^2{<p a \Pv\(p a ), 


(32) 
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which is akin to Eq. (31) except in the appearance of De rather than D in the upper limit of the summation. 
The third term (£3) concerns resonance. For the non-resonance condition, we could replace [16] 


e it( - b ' 3 E °-+ E -y E D = 6 a p5 a ^ + 6 a p6 70! - 5jp acT , 

from whence L 3 in Eq. (29) would reduce to 

L< f ] = E + E l I 1 ^/ 3 ) 1 2 + 'Y^{Va.\Pv\ ! Pc l ){Wl3\PAVp) 

OL 0^/3 

= 'y 1 I {^Pa \Pv Ity-’p) | T" 5 |-P; l^g) (V^/S \Py\V ft) ■ 

a=£0 a/3 

In the presence of resonance, however, other terms will also appear. 


(33) 


(34) 


L 3 = e lT< ' Efi Ea+E ' y E ° ) (Va.\PA L Pp)(Vo\PA l Pl) 


a.fi'ycr 


De 


(51 E e iT ^ G >) + 55 55 e it(g - G » ) )<^ a |P 1/ |^)(v3 a |P l/ |^ 7 ) - 55 KVc'lPvlVaW 

k: (a,/3)G^kl : m: (/3 ,ct)G Q m n: (a, 7 )G^ n a=l 


-De 


( 55 55 ^G k G, + 55 55 <5 G m G„)(<^a|-Pi,k/3)(^|P l -|^7) “ 55 K^I P ^l^a)| 2 (35) 

k: (a,/3)G^kl: (7,cr)6^l m: n: (a,7)Gt? n a=l 


£>e 

= 55(55 55 + 55 55 J^ai-P^x^i-P^) - 55 i<^«i p ^i^a)i 2 

k (a,/3)Gt/k (7,a)6^k (P,a-)Gt?k (a,7)e6k “ =1 

= X ( 3 nr) +55( 55 55 + 55 55 ^alP^X^alP^) 

k (a,@)eg k (7,c r )e^k|(a,^)^(7,cr) (/ 3 ,cr)e(a,7)6£ k |(P,c0^(a,7) 

= X 3 nr) +55 55 55 (^a|P«/|^)(^<T|P v |^7) 


(a,o’) GTk, 


f/3 7lef I (/3,7)#(a,a) 


= 4"'> + 4'> 


where 


■p 3 = E E E 


(^a | Pi/ 1 ^/ 3 ) (^o |P^ 1 ^ 7 ) > 


(36) 


(37) 


(a,<r)e5m s rr I (/3,7)^(“.< T ) 
(P,7)fe>m| (/3j7)7 4 (tTia) 


and in Eq. (35) we have used the identity 


e lT ( Gk Gl) = ^GkGi- 


(38) 


Noting that we allow resonance in the system, we do not necessarily have Sc k Gi — 'Ti- 
We now rewrite L [Eq. (30)] using Eqs. (32), (34), and (37) as 


D k 


’ D 


E 


\a=l 

|2 


l= E k^ip^i^)i 2 +(E(^i^i^)^§) + L 

<*#3=1 ' " 


(r) 

3 


dv 


< ;o max ^ i \(ip a \P u \<pi3)\ 2 + ^max^ i ( {ip a \Pv\<Pa) ~ ) + P 3 \ 


(39) 


(40) 


*7^/3(E{1,...,De} q:G{1,...,De} 

whose RHS is the very quantity X introduced in Eq. (26) in the strategy of the proof. We note that the expression for L here 

(r) 
' 3 


/ r \ 

differs from that in the non-degenerate-non-resonance case [16] in that here we have De rather than D, and a new term L\ 


has emerged. Equation (40) yields 


EpJ£]<E Pi/ max |(( / j q |P„|</? /3 )| 5 
q;//3g{1,...,De} 


- E p 


- q;G{1,...,De} 


max^ i ( (tp a \P v Wa) - -+ ) + E P „ [£+]. (41) 


■ (Oi 
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Following Ref. [16], the first two terms here can be bounded by using Eqs. (13) and (14) of Lemma 1. In order to set the scene 
to employ this lemma, we define the following orthonormal basis set for '-H\ 


l*«> 


\ip K )', 1 ^ K < D E 

14); d e + i^k^d, 


(42) 


where {|4}} are some vectors orthogonal to {|y> Q )}, chosen to complete the basis set, and thus (<E> K |<l> t ) = S KL . It is evident 
from Eq. (41) that 


E p v [X] sC Ep, 


max 




■ E p u 


max 

L«g{i,...,£>} 




\ _ d p 

D 


- Ep„ [X 


?] 


(43) 


Now we choose another orthonormal basis set for 9~( as {|£u l )}^. 1 such that P„ = \lo l )(oj l \, where % is the set of indices 

associated with the spectral representation of P v , with \%\ = d v , for all v £ {1 ,,M}. Hence if we define the unitary matrix 
U KL = ($ K |w t ), we can use Lemma 1 to obtain 


Ei 


^ ^ , 1(44 <pp)f 

La^/3G{1,...,-DE> 


■ E p 


max 
. q:G{1,...,-De} 


(444 - 




log D 9 d v log D 


D 


D 2 


Eq.p- 2 ) 10 log D 

< D ' 


(44) 


Now we show how to derive a bound on L\ and Ep, [L\ ] [in Eq. (41)]. From positivity of P„ (and hence P v eg) P v ), 

((4 <8 ({Per | - (dPfi | <8 <44 8 414 8 |4 ~ 14 8 |4) ^ 0, (45) 

we obtain 


Re [(<^c*| P i/ \(pfj)(ip a |Pi/|t^<y)] ^ „ ({ l fia\Pv\ i Pa} ( l Pa\Pv\ l fia) + {<P 0 \ Pv \ ^ ft) (^7 \Pv | <F 7 )) ■ 


(46) 


Combining Eqs. (37) and (46) yields 

■ (A .1 


4 r) ^E E E 


({ l Pa\Pv\ ( Pa} ((fieri Puller) T (v?/3 \Pv \ (<^7 | Pu | ^ 7 )) • (47) 


m (a,<j)€7 m -, p r I (/3.7)^(a,CT) 


For every pair (a, a) € > T m , the pair (8, 7 ) can have (f m — 2) different values. Furthermore, we have obtained two terms on 
the RHS of Eq. (47) which are equal. Thus 

L{ 3 < El(^ m “ 2 )] E (444(444 
<(D F - 2)E( E (444(444) 

m (a,rT)efrn 

= ( d f - 2)E(444(444 

olg 

= (D f - 2 ) E 444(444(444- (48) 

a a 

To get the third inequality we used the fact that each pair (a, a) appears only in one m-shell. It should be noted that the a = a 
case has automatically been included in our calculations, thus we do not need to consider it separately. Equation (48) gives 

Ep, [4 r) ] ^ (D f - 2) E P „ [E 4 2 4' 2 (444(444 

OLG 

= {D F - 2) E 444 E P u [(444(444]' 


OLG 


(49) 
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Since P v has a. fixed rank, averaging over random I\, can be equivalently replaced with averaging over random unitary operators 
U € SU (!)) according to the associated Haar measure, acting on a fixed projection P Uo 


E P „[4 r) ] =E U [L$\ P J 

< (D F -2)Y J \c a \ 2 \ca\ 2 ^u[^ a \U^P, 0 U\^ a ){^\U^P V0 U\i Pr7 )]. (50) 

OLG 


Since U is chosen based on the Haar measure of SU(-D), thus by applying random unitaries on | ip a ) (a specific/fixed vector in 
P,/-H) we indeed obtain vectors |£) which are uniformly-randomly distributed all over the D-dimensional Hilbert space 'Jd [36]. 
Hence we can rewrite Eq. (50) as 

■ 0 ) 


E P j4 r) ] < p F -2)^|c Q | 2 |c 7 | 2 E ? [||P l/0 |6)l| 2 ||^ol6) 


o: 7 


Eq. (9) 

< (D f -2)(E 5 [||PJ£)|| 4 ] +V 5 [||P, 0 |OH 2 ]) 


Lemma 2 


l«d„«-D 

< 


(d„\ 

/ + 1 \ 

\DJ 

VP + lJ 




It is evident that Ep^ [4^] becomes negligible if 


(D F -2)(^) 2 «L 


We insert Eqs. (44) and (51) into Eq. (41), which yields 

10 logD 


E P „[P]<^i^ + (P F -2)(^) 2 , 


(51) 


(52) 


(53) 


whose RHS is the very Ep y [X] introduced in the beginning of this subsection in the strategy of the proof. Next, following 
Eq. (27), we require that the RHS of Eq. (53) to be upperbounded by B5\ that is, 


or equivalently 


M 2 D r 10 log D 
S'E 2 dv - D 


M 2 r 10 log D 
SS's 2 i D 




< 6 , 




di/ 

< ~D 


(54) 


(55) 


On the other hand, in order to use Lemma 1 —in calculating the averages—we require condition (12). If we combine these 
conditions, we obtain 


max < C, 


10M 


2 r 


1 + 


(D f - 2 )dl 


SS's 2 l 10 D log D 


1 log D 
I D 


d v 1 

< 7) < C ' 


(56) 


V. SUMMARY 


It has been known that for the validity of the original quantum ergodic theorem, non-degeneracy and non-resonance properties 
for their Hamiltonians are among the conditions, which hold in typical systems. Here we, however, have considered systems 
which lack these conditions. We have proved a modified version of the quantum ergodic theorem which concerns validity of 
the ergodic hypothesis and normal typicality in atypical systems. We have shown that degeneracy does not considerably modify 
the condition of normal typicality or ergodicity, whereas the existence of resonance is more dominant for obstructing ergodicity. 
The effect of the non-resonance condition has been shown to come through degeneracy of energy sum structure of the system of 
question. 
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Appendix A: Proof that D\ = 2 for the non-resonant case 


The non-resonance condition implies that (see subsection IIB or Fig. 1) 


, 9 k = \Qk\ = 1, Vk ^ 0, 

whence Dq = 1. For example, if for some a , (3, 7 , and 0 we have 

(Al) 

E a — E 7 = Ep — E a = Gki 

for some k (or equivalently E a — Efj = Ey - E a = Gi), then 

(A2) 

(ct, 7 ) = (/3, ct) or (a, /3) = ( 7 , a), 

or equivalently 

(A3) 

<jk = {(a,7)}, 9k = 1, 

(A4) 

Q\ = {(a 5 /3)}, 91 = 1 - 

Equation (A2) can be recast in the following form too: 

(A5) 

Ea + E a = Ep + E y = F m , 

for some m, which means 

(A6) 

{(a, a), {a, a), (/3, 7 ), ( 7 ,^)} € J m . 

But the relations in Eq. (A3) both imply that 

(A7) 

(a,cr) = (/3,7). 

(A8) 


Note that there cannot be any other pair in f ln except those listed in Eq. (A7). Hence the non-resonance condition (Al) can be 
equivalently rewritten as 


/m=|Jm|=2, Vm, 


(A9) 


which in turn implies Dp = 2 . 

Remark .—Note that the existence of degeneracy contributes in our results in the paper mainly through modifying the estimates 
on the resulting resonance effects so that such effects become dominant. 


Appendix B: Calculating the average and variance of ||P„|yj Q ,}|| 2 with uniform measure on P„ 

Here we reproduce the average and variance of ||P„|<p Q ) || 2 , following the approach of Ref. [2]. Let us take {Iw,,)}^ as an 
orthonormal basis for the /9-dimensional Hilbert space .7/ such that 

= kX^I. (bd 

where % C {1,... ,D} constitutes the set of indices associated with P v , and J/„ = d v . We introduce another orthonormal 
basis {|$ K }fLi for 9 ~[related to {\tp a )Ya=\ as in Eq. (42). Now let us expand | ui L ) as 

D 

K) = + iy K )\$ K ), 

K,= 1 


(B2) 
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where the normalization condition (u> L |w t ) = 1 reads as ^^LiO 1 11 - 2 + Vk) = 1- This is the equation of the 2£)-dimensional unit 
hypersphere S 2D . In addition, we have 

\\PM\\ 2 = E(^ + ^)- (B3 > 

K,ey v 


Thus averaging over fixed-rank P„s with uniform distribution is equivalent to averaging uniformly over S 2D . This yields 


Ep„[||p,|^)II 2 ] = 2|X|e S2D [^], 


(B4) 


where E 52 D [Y] = \\/A 2 d\ fg 2 o dcr Y indicates the average over S 2D , in which dcr and A 2 D are the surface element and the 
total area of the unit 211-dimensional hypersphere. Straightforward calculations yield [2] 


E S 2D [ x k\ 


Vs 20 !® 2 ] 

Cs 222 [ x 2 k , x 2 ] 


1 

2D’ 

(D- 1 ) 
D 2 (D + iy 
1 

- D 2 {D + 1)’ 


K ^ i, 


(B5) 

(B 6 ) 

(B7) 


whence 


E pJ^k*)ll 2 ] =% 

\T I’ll E> I \ 112"i Eq-_(10) 1 

^pA\\p,M\\ } — 


fd u \ 2 (D — d u ) 
\D) (D + 1 ) ■ 


(B 8 ) 

(B9) 
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